Influence of diffusion on models for non-equilibrium wetting 
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It is shown that the critical properties of a recently studied model for non-equilibrium wetting 
are robust if one extends the dynamic rules by single-particle diffusion on terraces of the wetting 
layer. Examining the behavior at the critical point and along the phase transition line, we identify 
a special point in the phase diagram where detailed balance of the dynamical processes is partially 
broken. 
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I. INTRODUCTION 

In recent years there has been an considerable progress 
in the understanding of wetting transitions far from equi- 
librium. Such wetting processes are usually modelled 
by a non-equilibrium growth process of an interface in 
the vicinity of a hard-core wall which mimics an inert 
substrate. The dynamic rules of the growth process are 
chosen in such a way that detailed balance is violated, 
driving the system away from thermal equilibrium. De- 
pending on the balance between growth and evaporation 
rates the interface may either fluctuate in the vicinity 
or detach from the wall, leading to a wetting transition 
between a pinned and a moving phase. 

A simple model for non-equilibrium wetting was intro- 
duced in Ref. 1]. It is defined as a solid-on-solid growth 
process on a lattice with the restriction that the heights 
at neighboring sites may differ by at most one unit. While 
particles are deposited anywhere at the same rate, the 
rates for evaporation at the edges and in the middle of 
terraces are generally different. It was shown that such 
a growth process is effectively described by a Kardar- 
Parisi-Zhang equation 2] in a potential V: 
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Here £ is a white Gaussian noise, x denotes the position 
on the substrate, a controls surface tension, and V(h) is 
a potential of the form 



V(h) 
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voh if ft > 
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where v$ controls the propagation velocity of a free in- 
terface. The nonlinear term A[Vft(x, t)] 2 can be shown 
to break detailed balance and hence drives the system 
away from thermal equilibrium, even if the model is in 
the bound phase. 

As shown in Refs. 0, H, the properties of wetting 
transitions under non-equilibrium conditions differ signif- 
icantly from those at thermal equilibrium. For example, 
a non-equilibrium wetting transition is characterized by 
a different set of critical exponents. Moreover, if one in- 
troduces an additional short-range force between the sub- 
strate and the interface by adding a local well at ft = in 



the potential V(h), wetting models far from equilibrium 
may exhibit a phenomenon called phase coexistence, i.e., 
the bound and the moving phase are both thermodynam- 
ically stable in a certain region of the phase diagram. 

So far these theoretical predictions have not been ver- 
ified experimentally for various reasons. On the one 
hand it is very difficult to specify experimental conditions 
where detailed balance is broken on the level of micro- 
scopic transition rates. On the other hand, the models 
studied so far are unnatural in several respects. For ex- 
ample, in the model of Ref. [l[ the transition rates are 
assumed to be homogeneous while realistic experiments 
are always affected by inhomogeneities and impurities. 
Another idealistic assumption of the model, on which 
we will focus in the present work, is the circumstance 
that single atoms are not allowed to diffuse although 
they have a finite rate for evaporation. In experiments, 
however, one expects that the energy barrier for single- 
particle diffusion is of the same order or even less than the 
energy barrier for evaporation. Therefore, the question 
arises wether the phenomena of non-equilibrium wetting 
observed so far are robust against diffusion of solitary 
particles. 

Naively one may expect that the inclusion of single- 
particle diffusion leads to an effective term of the form 
DV 2 h in the continuum limit, which can be absorbed into 
the first term of the KPZ equation QJ without chang- 
ing the critical properties. However, in what follows we 
assume that only single particles are allowed to diffuse, 
while islands of two or more atoms remain stable. In 
this case the dynamics is changed in non-trivial way, in- 
cluding interesting phenomena such as partially broken 
detailed balance. 

The paper is organized as follows. In Sect. 2 we re- 
view the definition of the model introduced in 1] and 
modify the dynamical rules in order to incorporate single- 
particle diffusion. Sect. 3 summarizes essential proper- 
ties of equilibrium- and non-equilibrium wetting transi- 
tions. In Sect. 4 we first discuss a special case, where 
the dynamics of the model with diffusion is character- 
ized partially broken detailed balance, while we turn to 
the general behavior in Sect. 5. Finally our results are 
summarized in Sect. 6. 
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II. DEFINITION OF THE MODEL 

The original model introduced in Q is defined as fol- 
lows. The form of the interface is described by height 
variables hi = 0, 1, . . . at site i of a one-dimensional lat- 
tice with L sites, meaning that overhangs are not allowed. 
The configurations of the interface fulfill the so-called re- 
stricted solid-on-solid (RSOS) constraint, i.e., the height 
difference of two neighboring sites is restricted by 



\hi - hi 



±i 



< 1. 



(3) 



The model evolves random-sequentially according to the 
following dynamical rules. For each update with time 
increment At = ^, we randomly chose a site. Then one 
of the following processes is randomly selected according 
to specific probabilities (cf. Fig. QJ: 

a) An atom desorbs at the edge of a plateau with prob- 

ability jf\ 

hi -> min(hi-i, hi, h i+1 ) (4) 

b) A single atom desorbs with probability 

hi ^ hi - I (5) 

c) An atom desorbs from the interior of a plateau with 



probability jt\ 



hi — > hj 



1 



(6) 



d) An atom is deposited anywhere on the surface with 



probability jr\ 



hi^hi + l 



(7) 



Here N 



denotes the normalization constant. 



If the moves listed above would result into a configura- 
tion that violates the RSOS constraint Q or produces 
a height less than zero the move is abandoned, time is 
incremented, and the next site is selected. These rules 
can be generalized easily to higher dimensions. 

In this model there is no diffusion of solitary atoms 
along the surface. As outlined in the Introduction, one 
possible step to make the model more realistic is to in- 
troduce an additional rule for explicit diffusion. More 
specifically, a single atom is now allowed to diffuse a sin- 
gle step along a plateau at rate D by the following rule: 

e) Diffusion of a single atom on a plateau to the right or 



to the left with equal probability j^: 



hi — ► hi - 1 



hi 



±i 



hi 



±i 



1 



(8) 



Again this move is carried out only if the resulting con- 
figuration fulfills the RSOS constraint with non-negative 
heights. Note that this additional rule is not in conflict 
with the previous rules listed above. The normalization 
has to be replaced by N = D+2 l +p+q so that the proba- 
bilities change accordingly. 
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FIG. 1: Illustration of the dynamical rules. The transitions 
are: (a) desorption at an edge, (b) desorption of a single atom, 
and (c) desorption out of the middle of a plateau. Panel (d) 
displays three different situations where an atom can be de- 
posited. Furthermore the additional diffusive process consid- 
ered in the present work is illustrated in (e). 



III. PHENOMENOLOGY OF THE 
TRANSITION WITHOUT DIFFUSION 



Let us first review the properties of the wetting tran- 
sition in the original model without diffusion. Obviously, 
for a large growth rate q the interface propagates while 
for sufficiently small q it fluctuates near the wall. One 
can easily conclude that there is a critical growth rate 
q c where the propagation velocity of a free interface is 
zero. This defines a line of second-order phase transi- 
tions, which is shown in Fig.[2j This unbinding transition 
was shown to be continuous. 

Choosing p = 1 and q < 1 it was shown in [l| that 
the dynamical rules (a)-(d) obey detailed balance. This 
means that for any microscopic transition the reverse 
transition also exists and that the corresponding prob- 
ability currents in the stationary state compensate each 
other. Therefore, the bound interface approaches an 
equilibrium state in which the probability distribution 
of interface configurations is given by the Boltzmann en- 
semble. 

For p / 1, however, it can be shown that detailed 
balance is broken (see Fig. |2j). As before, the critical 
point q c can be determined by readjusting q in such a 
way that the asymptotic growth velocity is zero. 

A free interface is known to roughen algebraically in 
such a way that the width increases as 



w ~ t^f(t/L z 



(9) 



where / is a scaling function, L is the lateral system size, 
and (3 and z are the roughening and the dynamical expo- 
nent, respectively. It can be shown that for p = 1, where 
the bound phase exhibits detailed balance, the transition 
belongs to the so-called Edwards- Wilkinson universality 
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FIG. 2: Phase diagram of the model for non-equilibrium wet- 
ting introduced in Q in one dimension. A line of continuous 
phase transitions separates the pinned and the moving phase. 
For p = 1 the dynamical processes obey detailed balance (see 
text). Introducing diffusion shifts the phase transition line 
slightly but the overall structure of the phase diagram re- 
mains the same. 
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FIG. 3: Breaking of detailed balance for p ^ 1. The figure 
shows four possible configurations of the interface together 
with the corresponding transition rates. If Pa denoted the 
probability to find the interface in configuration A, detailed 
balance applied along the upper sequence of transitions would 
imply Pc = q 2 PA while the lower sequence of transitions 
would give Pc = (q 2 /p)Pa. Therefore, this counterexample 
demonstrates that detailed balance is broken for p ^ 1. 



exponent HHQ 

(1.18(1) forp<l 

3/4 for p = 1 . (10) 

0.23(1) forp>l 

This means that non-equilibrium wetting can be cate- 
gorized into three different universality classes, namely, 
the Edwards- Wilkinson class for p = 1 and two different 
KPZ-like classes depending on the sign of A in the KPZ 
equation, i.e., on p — 1. 

To each point of the phase transition line there is a 
profile P(h, t) that describes the probability to find the 
interface at height h. This probability distribution scales 
as 

P(h)=t-?g(ht-P), (11) 

where g is another universal scaling function describing 
the form of the profile. Therefore, profiles measured at 
different times can be collapsed onto one curve by plot- 
ting P(h)t@ against ht~^ . This method also allows one 
to estimate the exponent (3. As a major breakthrough, 
Prahofer and Spohn were able to compute the scaling 
function g for the case of a free interface analytically |?j . 
However, for non-equilibrium wetting (i.e. with a wall 
at zero height) the scaling function is still unknown al- 
though it could be approximated recently by mean-field 
techniques ^(J. In the following we investigate to what 
extent these results are robust with respect to diffusion 
of single particles. 



IV. TRANSITION WITH DIFFUSION FOR p = 1 

Let us now turn to the generalized model for non- 
equilibrium wetting which includes diffusion of single par- 
ticles (rule (e) in Fig. [I}. First we note that this process 
always breaks detailed balance, even for p = 1. This hap- 
pens when an atom diffuses towards the edge of an island 
where it sticks irreversibly. Obviously, this process can- 
not be reversed. So the interesting question is whether 
diffusion leads to a different type of critical behavior. 



class p , which is characterized by the exponents (3 = 1/4 
and z = 2. In the non-equilibrium case p / 1, how- 
ever, the asymptotic roughening behavior belongs to the 
KPZ universality class with the exponents (3 = 1/3 and 
z = 3/2 in one dimension Q. 

Introducing a wall at zero height induces an unbind- 
ing transition at q = q c , which may be characterized by 
the density of contact points po(t), where the interface 
has zero height. At the transition this quantity decays 
algebraically as po(t) ~ t _6> , giving rise to an additional 



-> 



FIG. 4: Breaking of detailed balance by single-particle dif- 
fusion. A particle diffuses to the edge of an island where it 
sticks irreversibly. Obviously, this process cannot be reversed 
so that detailed balance is broken. 

To answer this question we performed standard Monte 
Carlo simulations of a 1+1-dimensional system with L = 
4096 sites and T = 65536 time steps with and with- 
out diffusion, averaging over 100 independent runs. As 
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shown in Fig. [5J the width increases as t 1//4 in both 
cases. Surprisingly, detailed balance can be broken with- 
out changing the critical point q c = 1 and the exponent 
(3 = 1/4, i.e., without driving the transition away from 
the Edwards- Wilkinson universality class. However, as 
will be shown below, it is possible to understand this 
apparent contradiction. 
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FIG. 6: A possible configuration for a diffusive transition from 
site number 3 to site number 4. The arrows indicate those 
sites where atoms could adsorb or desorb, respectively. Obvi- 
ously, the diffusive move shown here increases the number of 
sites susceptible to deposition by 1. 
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FIG. 5: This graph shows that the width scales according to 
the Edwards- Wilkinson universality class as w ~ t 1//4 . 



FIG. 7: The mirror- symmetric configuration of the configu- 
ration shown in Fig. obtained by reflections h — > —h and 
x — > —x. Here a diffusive move increases the number of sites 
susceptible to desorption by 1, compensating the contribution 
of the configuration shown in Fig. El 



A. Invariance of the critical point for p = 1 

Let us first examine the possible transitions for each 
configuration of a free interface. For example, Fig. [S| 
shows a particular configuration of 6 sites where a parti- 
cle diffuses from site 3 to site 4. Initially we set D = 
and q = q c = 1, meaning that after sufficiently long time 
all allowed configurations occur with the same probabil- 
ity. Counting all sites where atoms can be deposited 
and similarly all sites from where particles can desorb we 
quantify how this segment of the interface contributes to 
the average growth velocity in the next instance of time. 
Since the asymptotic propagation velocity at the transi- 
tion is zero, these probabilities for deposition and evapo- 
ration averaged over all possible configurations have to be 
equal. Then, switching on diffusion, we pose the question 
whether diffusive moves destroy this balance or not. 

Indeed, for some configurations such as the one shown 
in Fig. El the number of potential adsorption and des- 
orption sites before and after the diffusive move changes. 
However, for each configuration there exists a mirror- 
symmetric configuration obtained by reflections h — » —h 
and/or x ^ x , as demonstrated in Fig. As can be 
seen, this mapping turns potential absorption sites into 
desorption sites and vice versa. Since both of them oc- 
cur with the same probability, the averaged change of 
the susceptibilities for adsorption and desorption cancel 
each other, meaning that the growth velocity remains 
zero even when diffusion is turned on. This argument 



proves that for p = 1 the transition point q c = 1 is inde- 
pendent of the diffusion rate D. 



B. Invariance of the height profiles for p — 1 

Taking a closer look at the violation of the detailed 
balance one can see that it is only a partial violation. 
First we note that diffusion does not change the total 
height H = hi, while all other transitions change H 
by one unit. Therefore, the configuration space of the 
system can be decomposed into sectors characterized by 
the value of H. Within those sectors diffusion is the only 
transition, while all other transitions take place between 
the sectors. Together with the previous finding that the 
total rates for desorption and adsorption between the sec- 
tors is invariant we conclude that the probability to find 
the system in a particular sector does not depend D. 
Moreover, it is obvious that diffusive moves within the 
sectors do not change the height profile P(h). This means 
that the height profile does not depend on the value of 
D and the dynamics between the sector still obeys de- 
tailed balance. Therefore, the present model provides an 
interesting example of a system where detailed balance 
is partially broken. 
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V. TRANSITION WITH DIFFUSION FOR p / 1 

For p 7^ 1 we expect that the critical point q c depends 
on the diffusion rate. For example, for the configurations 
shown in Fig. |S| and Fig. {7\ the aforementioned compen- 
sation mechanism is no longer possible since rates for 
deposition q < 1 and desorption at the edges r — 1 are 
different. In fact, considering all possible configurations 
of a segment of six sites, one finds that a diffusive move 
creates 7 additional sites where an atom could desorb 
with rate p and 6 additional sites where an atom could 
be deposited with rate q while there will be one site less 
from where an atom could desorb with rate r. As the 
rates in the nonequilibrium case along the phase transi- 
tion line are not equal any more the phasetransition line 
with diffusion is different from the one without diffusion. 




(1/3) 



h/t 



FIG. 8: Data collapse for p — 0.3 for different times 



In order to verify whether diffusion changes the critical 
behavior we determined interface profiles at the critical 
point. As mentioned in Sect. Hill different curves of pro- 
files recorded at different times can be collapsed onto a 
single curve by using the scaling form (|TT|) with the cor- 
rect scaling exponent (3. For D = 1 we find that the 
critical point does not change significantly. Collapsing 
our profiles for with the same value of (3 we obtain a con- 
vincing data collapse, supporting that for p ^ 1 diffusion 
does not change the universality class of the transition. 



1) Single particles on a plateau continue to diffuse at 
infinite rate (unable to leave the island) until an- 
other atom is deposited on the plateau. This diffu- 
sion process is not simulated explicitly, instead the 
nucleation site can be chosen randomly in agree- 
ment with the definition of the model. 

2) Single particles diffusing on a terrace instanta- 
neously attach at the neighboring island. 

3) Single particles in a well will instantaneously stick 
randomly to one of the neighboring islands with the 
probabilities 



Pi 



d r 



di + d r 



Pr 



di 



di + d r 



(12) 



where d\ and d r denotes the distances to the left 
and the right island, respectively. Note that this 
linear law holds only in the one-dimensional case. 

Looking at typical snapshots of the interface for D = 
and D = oo in Fig.[2|there is no big difference apart from 
the fact that in the latter case there are no single atoms 
any more. Carrying out numerical simulations similar to 
those presented in the previous sections we obtain similar 
results, in particular, the height profiles for D = and 
D — > oo can be collapsed onto a single curve in each of 
the three cases p < 1, p = 1, and p > 1. This means 
that the critical exponents do not depend on D so that 
the universality class of the transition is not changed 
by single-particle diffusion, even in the limit D — > oo. 
This result is plausible since the arguments discussed in 
Sects. HV AIIV Bl do not depend on the value of D. 




site i 



VI. BEHAVIOR IN THE LIMIT D -> oo 

So far we discussed the behavior of the extended model 
for finite values of the diffusion rate D. How robust are 
the results when taking D to infinity? To study this 
limit efficiently, we reformulate the dynamic rules with- 
out changing the definition of the model as follows. After 
each update we identify isolated particles that could pos- 
sibly diffuse. Each of them is found in one of the following 
local configurations and evolves accordingly: 



FIG. 9: Snapshots of a free interface after 100 time steps 
with infinite and without diffusion. The profile with diffusion 
is characterized by the absence of single atoms. 



VII. CONCLUSIONS 

In this paper we have generalized a previously intro- 
duced model for non-equilibrium wetting in such a way 
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that freshly deposited adatoms are allowed to diffuse. 
This generalization is necessary to compare this model 
with possible experimental applications, where the en- 
ergy barrier for horizontal diffusion is usually of the same 
order or even less than the energy barrier for evapora- 
tion. As a main result we observe that single-particle 
diffusion in 1+1 dimensions does not change the critical 
behavior at the transition, i.e., all predictions made for 
non-equilibrium in previous works remain valid. 

The results for the special point p = 1, where the origi- 
nal model exhibited detailed balance in the bound phase, 
are particular surprising. Although diffusion breaks de- 
tailed balance and even changes the visual appearance of 
the interface (see Fig. the transition point q c = 1 as 
well as the rescaled height profile remain invariant, even 
in the limit D — » oo. This apparent contradiction can be 
resolved by observing that detailed balance is only par- 



tially broken within sectors of constant integrated height. 
This example demonstrates that broken detailed balance 
does not always drive the system away from equilibrium. 

Extending models for non-equilibrium wetting by 
single-particle diffusion is only a first step towards a more 
realistic modelling aiming at a quantitative comparison 
with possible experiments. In addition, a detailed survey 
of energy barriers and the corresponding transition rates 
under typical experimental conditions would be needed. 
Moreover, it would be interesting to examine the influ- 
ence of disorder caused by inhomogeneities and defects 
of the substrate. 
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